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Abstract 

Let p > 3 be a prime. Euler numbers E p s first appeared in H. S. Vandiver's work 
(1940) in connection with the first case of Fermat Last Theorem. Vandiver proved that 
x p _|_ yP — z p has no solution for integers x,y,z with gcd(xy z , p) — 1 if E p _s = 
(modp). Numerous combinatorial congruences recently obtained by Z.-W. Sun and 
by Z.-H. Sun involve the Euler numbers E p ^. This gives a new significance to the 
primes p for which E p s = (modp). 

For the computation of residues of Euler numbers B p _3 modulo a prime p, we use 
the congruence which runs significantly faster than other known congruences involv- 
ing E p -3. Applying this congruence, a computation via Mathemat ica 8 shows that 
only three primes less than 10 7 satisfy the condition E p s = (modp) (such primes 
are 149, 241 and 2946901, and they are given as a Sloane's sequence A198245). By 
using related computational results and statistical considerations similar to those on 
search for Wieferich and Fibonacci- Wieferich and Wolstenholme primes, we conjec- 
ture that there are infinitely many primes p such that E p s = (modp). Moreover, 
we propose a conjecture on the asymptotic estimate of number of primes p in an inter- 
val [x, y] such that S p _3 = A (modp) for some integer A with \A\ G [K, L]. 

Keywords: Euler number, E p s, congruence modulo a prime, supercongruence, 
Fermat quotient 

1. Introduction 

Euler numbers E n (n — 0,1,2,...) (e.g., see 1 1 3. pp. 202-203]) are integers 
defined recursively by 



E = 1, and y\E n ^ k 

0<k<n \ ' 



5 for n = l,2,3,... 
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(it is well known that EWi-i = for each n = 1,2,...)- The first few Euler numbers 
are E a = l,E 2 = -l,E 4 = 5,£ 6 = -61, E* = 1385, E w = -50521, E 12 = 
2702765, E u = -199360981, E 16 = 19391512145. It is well known that Euler num- 
bers can also be defined by the generating function 



n=0 



It is well known that E n = E n (0) (n — 0,1,.. .), where E n (x) is the classical Euler 
polynomial (see e.g., 61 et seq.]). 

Recall that Bernoulli numbers B n (n — 0,1,2,.. .) are rational numbers defined 
by the formal identity 



n=0 



It is easy to see that B n = for odd n > 3, and the first few nonzero terms of (B n ) 
are B = 1, B x = -1/2, B 2 = 1/6, B 4 = -1/30, B 6 = 1/42 and B$ = -1/30. It is 
well known that B n = B n (Q), where B n (x) is the classical Bernoulli polynomial (see 



e.g., 11151, p. 61 et seq.]). 

A significance of Euler numbers, and especially of -B p _3 with a prime p, is closely 
related to Fermat Last Theorem (see Lecture X, Section 2]). In 1850 Kummer 
(see e.g., iQJ, Theorem (3A), p. 86 and Theorems (2A)-(2F), pp. 99-103] proved that 
Fermat Last Theorem holds for each regular prime, that is, for each prime p that does 
not divide the numerator of any Bernoulli number B 2n with n = 1, 2, . . . , (p — 3)/2. In 
1940 H. S. Vandiver [24] likewise proved for Euler-regular primes. Paralleling the pre- 
vious definition of a (irr)regular prime (with respect to the Bernoulli numbers) follow- 
ing Vandiver M24I1 . a prime p is said to be Euler-irregular primes (shortly E-irregular) 
if and only if it divides at least one of the Euler numbers E2 n with 1 < n < (p — 3)/2. 
Otherwise, that is if p does not divide E%, E4, . . . , E p _3, a prime p is called E-regular. 
The smallest £ , -irregular prime is p = 19, which divides ^(lO) = —50521. The first 
few E-irregular primes are 19, 31, 43, 47, 61, 67, 71, 79, 101, 137, 139, 149, 193, 223, 241 
(with p = 241 dividing both E210 and an d hence having an E'-irregularity index 
of 2) (see [J]). In 1954 L. Carlitz fljj] proved that there are infinitely many E'-irregular 
primes p, i.e., p \ E2E4 ■ ■ ■ E p ^. Using modular arithmetic to determine divisibility 
properties of the corresponding Euler numbers, the i?-irregular primes less than 10000 
were found in 1978 by R. Ernvall and T. Metsiinkyla 0]. 



In his book fl 1 3L p. 203] P. Ribenboim noticed that "it is not all surprizing that 
the connection, via Rummer's theorem, between the primes dividing certain Bernoulli 
numbers and the truth of Fermat's theorem, would suggest a similar theorem using the 
Euler numbers." Vandiver [24] proved that x p + y p — z p has no solution for integers 
x, y, z with gcd(xyz,p) = 1 if -Ep-3 = (modp). The analogous result was proved 
by Cauchy (1847) and Genocchi (1852) (see £[1 p. 29, Lecture II, Section 2]) with the 
Bernoulli number -Bp-3 instead of -E p _3. Further, in 1950 M. Gut [8] proved that the 
condition E p _s = E p _§ = E p ^j = E p ^g = E p -\\ = (modp) is necessary for the 
Diophantine equation x 2p + y 2p = z 2p to be solvable. 

Furthermore, numerous combinatorial congruences recently obtained by Z.-W. Sun 



in I201- II22I1 and by Z.-H. Sun in [17] involve Euler numbers E p -z with a prime 
p. Many of these congruences become "supercongruences" if and only if E p s = 
(modp) (A supercongruence is a congruence whose modulus is a prime power.) 
This gives a significance to primes p for which E p ^ = (modp). The first two 
primes 149 and 241 have also been discoverded by Z.-W. Sun 112011 . 
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In this note, we focus our attention to the computational search for residues of Euler 
numbers S p _3 modulo a prime p. By the congruence obtained in 1938 by E. Lehmer 
p. 359], for each prime p > 5 

b/4] 

J2 tj = (-l) (p - 1)/2 4£; p - 3 (modp), (1) 
fc=i 

where [a] denotes the integer part of a real number a. Usually (cf. fl), if £„_ 3 = 
( mod p) then we say that (p, p — 3) is an E-irregular pair. It was founded in 14J that 
in the range p < 10 4 (p,p — 3) is an _E-irregular pair for p = 149 and p = 241. 

For our computations presented in Section 3 we do not use Lehmer's congruence (1) 
including harmonic number of the second order. Our computation via Mathemat i ca 
8 which uses the expression including the harmonic number (of the first order) is 
very much faster than those related to the congruence (1). Here we report that only 
three primes less than 10 7 satisfy the condition E p -3 = (modp). Using our com- 
putational results and statistical considerations similar to those in relation to a search 
for Wieferich and Fibonacci-Wieferich and Wolstenholme primes (cf. |2, p. 447] 
and iUHd . we conjecture that there are infinitely many primes p such that E p _3 
(modp). 



2. A congruence used in our computation 

Here, as usually in the sequel, for integers m, n, rs with « / and s ^ 0, and a 
prime power p a we put m/n = r/s (mod p e ) if and only if ras = nr (mod p e ), and 
the residue class of m/n is the residue class of ran' where n' is the inverse of n modulo 
P e - 

In what follows p always denotes a prime. The Fermat Little Theorem states that if 
p is a prime and a is an integer not divisible by p, then a p_1 = 1 (mod p). This gives 
rise to the definition of the Fermat quotient ofp to base a, 

qP- 1 - 1 
QpW : = 1 

which is an integer. It is well known that divisibility of Fermat quotient q p (a) by p 
has numerous applications which include the Fermat Last Theorem and squarefreeness 
testing (see 01 > 01 an d iH)- If <7p(2) is divisible by p, p is said to be Wieferich prime. 
Despite several intensive searches, only two Wieferich primes are known: p = 1093 
and p = 3511 (see |0] and fl). Another class of primes initially defined because 
of Fermat Last Theorem are Fibonacci-Wieferich primes, sometimes called Wall-Sun- 
Sun primes. A prime p is said to be Fibonacci-Wieferich prime if the Fibonacci number 
^ p zii) * S ^' v ^ s ^^ e ky P 2 ' wnere (f) denotes the Legendre symbol (see 11 810 . A search 
in iTTlfl and shows that there are no Fibonacci-Wieferich primes less than 9.7 x 10 14 . 

For the computation of residues of Euler numbers £^-3 modulo a prime p, it is 
suitable to use the following congruence which runs significantly faster than Lehmer's 
congruence (1). 

Theorem (Q3, Theorem 4.1 (iii)]). Let p > 5 be a prime. Then 

[P/4] 

£ i + 3q p (2) - fq p (2 f ee (-l^/V^-S (mod p\ (2) 
k=l 

where [a] denotes the integer part of a real number a. 



3 



Proof. Quite recently, Z.-W. Sun [20, Proof of Theorem 1.1, the congruence after (2.3)] 
noticed that by a result of Z.-H. Sun [ 17, Corollary 3.3], 

^^ P ( 2 )^M2) 2 -(-l) (?+1)/ H-3 (mod/). (3) 

fc=i 

On the other hand, we have 



E ^— E I- 2 E i- E I- E ■ 

fe=l fc=l l<Jo<(p-l)/2 fe = l j=l J 



By the classical congruence proved in 1938 by E. Lehmer [9J the congruence (45), p. 
358], for each prime p > 5 

(p-1)/2 1 

J2 ^ = -2q p (2)+pq p (2) 2 (mod/). (5) 
fc=l 

Substituting the congruence (5) into (4), we obtain 



(p-1)/2 f _ llfe _i j [p/4] 

E M = -2q p (2)+pq p (2) 2 (modp 2 )--J2- (mod p 2 ). (6) 

k=l 3=1 J 



Finally, substituting (6) into (3), we immediately obtain (2). □ 



3. The computation 

Using the congruence (2), a computation via Mathematica 8 shows that only 
three primes less than 10 7 satisfy the condition -E p _3 = (modp) (such primes are 
149, 241 and 2946901, and they are given as a sequence A198245 in H). Notice 
also that in 201 1 11121 p. 3, Remarks], the author of this article reported that these three 
primes are only primes less than 3 x 10 6 . 

Recall that investigations of such primes have been recently suggested by Z.-W. Sun 
in |2(J; namely, in luOi Remark 1.1] Sun found the first and the second such primes, 
149 and 241, and used them to discover curious supercongruences (1 .2)— (1 .5) from 
Theorem 1.1 in [20] involving E p s. 

Motivated by search for Wieferich and Fibonacci-Wieferich primes given in [2[] and 
J3| and search for Wolstenholme primes given in [11], here we use similar computa- 



tional considerations for Euler numbers -Ep_3 where p is a prime. Our computational 
results presented below suggest two conjectures on numbers i? p _3 that are analogous to 
those on Wieferich (101, J3D) and Wolstenholme primes lUUl . Accordingly, we search 
primes p in the range [10 5 ,5 x 10 6 ] such that E p ^ 3 = A(modp) with \A\ < 100 
and/or 10 4 • \A/p\ < 1. Our search employed the congruence (2) which runs signifi- 
cantly faster than Lehmer's congruence (1) and than the code 

Print [{Prime [n] }, Mod [EulerE [Prime [n] -3 ] , Prime [n] ] ] 

Here EulerE [k] gives Ef. and Mod [ a, m] gives a(mod m). 
in Mathemat ica 8, as well as than some other known congruences involving Euler 
number E p ^ 3 . 
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Namely, in order to obtain data of Table 1 below concerning primes p with 10 5 < 
p < 5 x 10 6 we used the code: 

Do [ If [Max [Min [Mod [ (Mod [Numerator [HarmonicNumber [Floor [Prime [n] /4] ] ] , 
Prime [n] " 2 ] PowerMod [Denominator [HarmonicNumber [Floor [Prime [n] /4 ] ] ] , 
-1, Prime [n] "2 + 3* (2 * (Prime [n] -1) -1) /Prime [n] -PowerMod [2, -1, Prime [n] "2] 
* (3*Prime [n] ) * ( (2" (Prime [n] -1) /Prime [n] *2) / ( (-1) * ( (Prime [n] +1) /2) 
*Prime [n ] ) , Prime [n] ] , Prime [n] -Mod [ (Mod [Numerator [HarmonicNumber 
[Floor [Prime [n] / 4 ] ] ] , Prime [n] " ] *PowerMod [Denominator [HarmonicNumber 
[Floor [Prime [n] /4 ] ] ] , -1, Prime [n] "2 + 3* (2 " (Prime [n] -1) -1) /Prime [n] 
-PowerMod [2, -1, Prime [n] "2 ] * (3*Prime [n] ) * ( (2 " (Prime [n] -1) -1) /Prime [n] ) "2) 
/ ( (-1) " ( (Prime [n] +1) /2) *Prime [n] ) , Prime [n] ] ] ] ==1000, Print [{n}, 
{Prime [n] }, {Mod [ (Mod [Numerator [HarmonicNumber [Floor [Prime [n] /4] ] ] , 
Prime [n ] " 2 ] * PowerMod [Denominator [HarmonicNumber [Floor [Prime [n]/4]]],-l. 
Prime [n] "2 + 3* (2" (Prime [n] -1) -1) /Prime [n] -PowerMod [2, -1, Prime [n] "2] 
* (3*Prime [n] ) * ( (2" (Prime [n] -1) /Prime [n]"2)/((-l)"( (Prime [n] +1) /2) 
*Prime [n] ) , Prime [n] ] }, {Prime [n] -Mod [ (Mod [Numerator [HarmonicNumber 
[Floor [Prime [n] /4 ] ] ] , Prime [n] "2 ] *PowerMod [Denominator [HarmonicNumber 
[Floor [Prime [n] /4 ] ] ] , -1, Prime [n] "2+3* (2" (Prime [n] -1) -1) /Prime [n] 
-PowerMod [2, -1, Prime [n] "2 ] * (3*Prime [n] ) * ( (2" (Prime [n] -1) /Prime [n] "2) 
/ ( (-1) ~ ( (Prime [n] +1) 12) *Prime [n] ) , Prime [n] ] }] ] , {n, i, j}] 

Here Mod [ a, m] gives a(modm), PowerMod [a, b, m] gives a b (mod m) (and 
is faster than Mod [a"b, m] . 

Further, in order to verify that there are no primes p between 5 x 10 6 and 10 7 such 
that -Ep-3 = (modp), we used the following code which is very much faster the 
previous code: 

Do [Print [{n}, {Prime [n] }, Mod [Numerator [2*HarmonicNumber [Floor [ 
Prime [n]/4]]+6*(2" (Prime [n] -1) -1) /Prime [n] -3* (2" (Prime [n] -1) -1) "2 
/Prime [n] ] , Prime [n]"2]],{n,i,j}] 

Certainly A = A(p) can take any of p values (modp). Assuming that A takes 
these values these values randomly, the "probability" that A takes any particular value 
(say 0) is 1/p. From this, in accordance to the heuristic given in H related to the 
Wieferich primes, we might argue that the number of primes p in an interval [x, y] such 
that £" p -3 = (modp) is expected to be 



If this is the case, we would be only expect to find about 0.998529(~ 1), such primes 
in the interval [10 7 , 10 19 ]. On the other hand, since 9999991 is the greatest prime 
less than 10 7 and is is actually 664589th prime, by above estimate, we find that in 
the interval [2, 10 7 ] we can expect about X^2<p<io 7 VP = Y^ktli 89 ^/Pk ~ 3.04145 
primes p such that E p s = (modp) (pk is a fcth prime); as noticed previously, our 
computation shows that all these primes are 149, 241 and 2946901. 




(7) 
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Table 1. Primes p with 10 5 < p < 5 x 10 6 for which E p _3 = A (modp) with 
\A\ < 100 and/or with related values \ A/p\ < 10 -4 (given in multiples of 10~ 4 ) 



V 


A 


\A/p\ 


V 


A 


\A/p\ 


105829 


-74 


> 1 


1355269 


-60 


0.442717 


111733 


45 


> 1 


1392323 


-29 


0.208285 


127487 


38 


> 1 


1462421 


-78 


0.533362 


130489 


-27 


> 1 


1546967 


-43 


0.277963 


131617 


9 


0.683802 


1743271 


107(> 100) 


0.613789 


162847 


-85 


> 1 


1794049 


-131(< -100) 


0.730192 


165157 


-46 


> 1 


1808497 


-121(< -100) 


0.669109 


171091 


- 17 


0.993623 


1952131 


-153(< -100) 


0.783759 


171449 


7 


0.408285 


1986539 


-157(< -100) 


0.790319 


191237 


37 


> 1 


2053873 


18 


0.087639 


192961 


63 


> 1 


2114251 


211(> 100) 


0.997989 


200461 


7 


0.349195 


2236349 


4 


0.017886 


209393 


27 


> 1 


2342381 


143(> 100) 


0.610490 


245471 


39 


> 1 


2410627 


-219(< -100) 


0.908477 


246899 


-54 


> 1 


2472731 


230(> 100) 


0.930146 


276371 


-69 


> 1 


2583011 


159(> 100) 


0.615561 


290347 


10 


0.344415 


2619847 


224(> 100) 


0.855011 


292183 


53 


> 1 


2740421 


225 


0.821042 


306739 


-42 


> 1 


2890127 


-34 


0.117642 


317263 


-35 


> 1 


2946901 








321509 


84 


> 1 


3279833 


-111(< -100) 


0.338432 


342569 


25 


0.729780 


3290689 


200(> 100) 


0.607775 


422789 


-40 


0.946098 


3312653 


228(> 100) 


0.688270 


429397 


-62 


> 1 


3340277 


226(> 100) 


0.676591 


440047 


82 


> 1 


3355813 


116(> 100) 


0.345669 


479561 


31 


0.646425 


3652613 


-290(< -100) 


0.793952 


501317 


60 


> 1 


3818131 


-318(< -100) 


0.832868 


546631 


92 


> 1 


3852677 


75 


0.194670 


628301 


73 


> 1 


3960377 


-48 


0.121201 


636137 


25 


0.392997 


4007747 


190(> 100) 


0.474082 


656147 


-68 


> 1 


4121503 


-270(< -100) 


0.655101 


659171 


-22 


0.333753 


4171229 


153(> 100) 


0.366798 


687403 


-4 


0.058190 


4343659 


-252(< -100) 


0.580156 


717667 


-42 


0.585230 


4392007 


55 


0.125227 


719947 


53 


0.736165 


4418497 


70 


0.158425 


766261 


-8 


0.104403 


4475707 


193(> 100) 


0.431217 


801709 


53 


0.661088 


4541501 


120(> 100) 


0.264230 


920921 


-82 


0.890413 


4551973 


-362(-100) 


0.795260 


924727 


-8 


0.086512 


4564939 


-63 


0.138008 


1064477 


106(> 100) 


0.995794 


4631399 


367(> 100) 


0.792417 


1080091 


42 


0.388856 


4674347 


302(> 100) 


0.646080 


1159339 


-38 


0.327773 


4706047 


220(> 100) 


0.467484 


1202843 


21 


0.174586 


4751599 


-279(< -100) 


0.587171 


1228691 


15 


0.122081 


4761677 


200(> 100) 


0.420020 


1285301 


47 


0.365673 


4869517 


-100 


0.205359 


1336469 


-5 


0.037412 


4898099 


-236(< -100) 


0.481820 


1353281 


78 


0.576377 


4928503 


-173(< -100) 


0.351019 
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The second column of Table 1 shows that there are 61 primes between 10 5 and 
5 x 10 6 for which \A\ < 100. Since the "probability" that \A\ < 100 for a prime 
p 3> 200 is equal to 201/p, it follows that expected number of such primes between 
Mth prime p M and Nth prime p N with N > M > 1000 (that is, p N > p M > 1000) 
is equal to 

Q(N, M, 100) = 201 W 

PM<P<PN P 

where the summation ranges over all primes p such that pu < P < Pn- In particular, 
for the values M = 9593 and N = 348513 which correspond to the interval [10 5 , 5 x 
10 6 ] containing all primes from Table 1, we have 

Q(348513, 9593, 100) = 201 ^ - « 201 • 0.292251 = 58.742451. (9) 

10 5 <p<5xl0 6 P 

On the other hand, Table 1 shows that there are 61 primes between 10 5 and 5 x 10 6 
for which \A\ < 100, which is « 3.8431% greater than related "expected number" 
58.742451. 

Because our program recorder all p with "small \ A\", that is, with \ A\ < 100, we 
compiled a large data set which can be used to give more rigorous (experimental) con- 
firmation of both our Conjectures 1 and 2. Indeed, our program recorded 568 primes p 
in the interval [10 5 , 5 x 10 6 ] for which \ A\ < 1000. On the other hand, according to 
the formula (9), it follows that expected number of such primes is equal to 

Q(348513, 9593, 1000) = 2001 ^ ~ ~ 2001 ' 0.292251 = 584.794251 

10 5 <p<5xl0 6 P 

(10) 

which is w 2.956% greater than related "expected number" 568. 

Instead, of selecting values based on |^4| < 100, we suggest to select them based 
on A/p < q x 10~ 4 (e.g., q = 1) that would be consistent with the original selection 
criterion. In particular, in the third column of Table 1 there are 72 primes p contained 
in the interval [10 5 , 5 x 10 6 with related values 10 4 x A/p < 1. 

Furthermore, since the "probability" that \ A/p\ < 10~ 4 for a prime p ^> 10000 is 
equal to 

2 [loooo] + 1 _ 2 
p ~ 10000' 

it follows that expected number of such primes between Mth prime pu and iVth prime 
p N with N > M > 1000 (that is, p N > p M > 10000) is equal to 

2(N - M) 



P{N, M) = 



10000 



In particular, for the values M = 9593 and N = 348513 which correspond to the 
range (10 5 , 5 x 10 6 ) of all primes from Table 1, we have 

x 677840 
P(348513,9593) = — — = 67.7840 
v ' ' 10000 

which is w 5.855% less than 72. 

All the previous considerations and the well known fact that the series 



E 1 

p prime 
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diverges suggest the following conjecture. 

Conjecture 1. There are infinitely many primes p such that E p _s = (modp). 
Since 

V - w log log x - log log y, 

in view of the previous comparison of our computational results with expected number 
of primes p E [10 5 , 5 x 10 6 ] for which < 100 given by (9) (or primes p £ 

[10 5 , 5 x 10 6 ] for which < 1000 given by (10)), we can assume that expected 

number of primes p in an interval [x, y] such that K < \A(p)\ < L is asymptotically 
equal to (cf. (7)) 

2(L- K)- (log log b -log log a). (11) 

Using a larger data set which our program recorded, consisting of total 568 pairs 
(p,A(p)) such that p e [10 5 ,5 x 10 6 ] and \A(p)\ < 1000, we obtain experimental 
results presented in Table 2. In Table 2 the values in "column fc" and in first and second 
row reflect the number of p 6 [10 5 , 10 6 ] and p 6 [10 6 , 5 x 10 6 ], respectively, such that 
A = A(p) e [k x 100, (k + 1) X 100] (k = 0, 1, . . . 9). Expected numbers given in the 
last column of Table 2 are calculated by the formula (11). 



Table 2. 





k 




Interval 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Expected 


[10 b ,10 B ] 


42 


51 


37 


30 


29 


24 


31 


34 


42 


44 


36.464 


[10 6 ,5 x 10 6 ] 


22 


23 


26 


20 


22 


22 


21 


24 


21 


20 


22.039 



Table 2 presents a small snapshot of our experimental results. Notice that by the 
data of the last row, the relative error between the conjectured and experimental values 
for k = 0, 1, . . . , 9 are respectively equal to 0.18%, 4.18%, 15.23%, 10.20%, 0.18%, 
0.18%, 4.95%, 8.17%, 4.95%, 10.20%. Accordingly, we propose the following con- 
jecture (cf. the same conjecture in [ 3, Conjecture 6.1] concerning the Wieferich primes; 
see also J2, Section 3]). 

Conjecture 2. The number of primes p G [a,b] such that \A\ = \A(p)\ G [K, L] is 
asymptotically 

2 (L — K) ■ (log log b — log log a) . 

Remarks. Recall that a prime p is said to be a Wolstenholme prime if it satisfies the 
congruence 

(2p - 1 
VP — l 

or equivalently (cf. ifioi Corollary on page 386]; also see fl) that p divides the numera- 
tor of i? p -3. The only two known such primes are 16843 and 2124679, and by a result 
of R.J. Mcintosh and E.L. Roettger from 11 ll pp. 2092-2093], these primes are the 
only two Wolstenholme primes less than 10 9 . Nevertheless, using similar arguments to 
those given in Section 3 of this paper, Mcintosh ifTot page 387] conjectured that there 
are infinitely many Wolstenholme primes. 



J = 1 (mod p 4 ), 



8 



References 



References 

[1] L. Carlitz, Note on irregular primes, Proc. Amer. Math. Soc. 5 (1954) 329-331. 

[2] R. Crandall, K. Dilcher and C. Pomerance, A search for Wieferich and Wilson 
primes, Math. Comp. 66 (1997) 443-449. 

[3] KG. Dorais and D. Klyve, A Wieferich prime search up to 6.7 x 10 15 , J. Integer 
Seq. 14(2011) Article 11.9.2. 

[4] R. Ernvall and T. Metsankyla, Cyclotomic invariants and i?-irregular primes, 
Math. Comp. 32 (1978) 617-629. 

[5] R. Ernvall and T. Metsankyla, On the p-divisibylity of Fermat quotients, Math. 
Comp. 66(1997) 1353-1365. 

[6] J. W.L. Glaisher, On the residues of the sums of products of the first p— 1 numbers, 
and their powers, to modulus p 2 or p 3 , Q. J. Math. 31 (1900) 321-353. 

[7] A. Granville, Some conjectures related to Fermat' s Last Theorem, Number Theory 
(Banff, AB, 1988), de Gruyter, Berlin, 1990, 177-192. 

[8] M. Gut, Eulersche Zahlen und grosser Fermat'scher Satz, Comment. Math. Helv. 
24 (1950) 73-99. 

[9] E. Lehmer, On congruences involving Bernoulli numbers and the quotients of 
Fermat and Wilson, Ann. Math. 39 (1938) 350-360. 

[10] R.J. Mcintosh, On the converse of Wolstenholme's theorem, Acta Arith. 71 (1995) 
381-389. 

[11] R.J. Mcintosh and E.L. Roettger, A search for Fibonacci-Wieferich and Wolsten- 
holme primes, Math. Comp. 76 (2007) 2087-2094. 

[12] R. Mestrovic, An exstension of a congruence by Kohnen, 13 pages, preprint 
larXiv: 1109 . 2340t v/3 [math. NT] (2011). 

[13] P. Ribenboim, 13 Lectures on Fermat's Last Theorem, Springer- Verlag, New 
York, Heidelberg, Berlin, 1979. 

[14] N.J. A. Sloane, Sequence A198245 in OEIS (On-Line Encyclopedia of Integer 

Sequences), |http : //oeis . org/A1982 4 5| 

[15] H.M. Srivastava and J. Choi, Series Associated with the Zeta and Related Func- 
tions, Kluwer Academic Publishers, Dordrecht, Boston and London, 2001. 

[16] Z.-H. Sun, Congruences concerning Bernoulli numbers and Bernoulli polynomi- 
als, Discrete Appl. Math. 105 (2000) 193-223. 

[17] Z.-H. Sun, Congruences involving Bernoulli and Euler numbers, J. Number The- 
ory, 128 (2008) 280-312. 

[18] Z.-H. Sun and Z.-W. Sun, Fibonacci numbers and Fermat's last theorem, Acta 
Arith. 60(1992) 371-388. 



9 



[19] Z.-W. Sun, Binomial coefficients, Catalan numbers and Lucas quotients, 
Sci. China Math. 53 (2010) 2473-2488; preprint larXiv: 90 9 . 5 64 81 ^11 
[math. NT] (2010). 

[20] Z.-W. Sun, On Delannoy numbers and Schroder numbers, J. Number Theory 131 
(201 1) 2387-2397; preprint larXiv :1009.2486V 4 [math. NT] (2011). 

[21] Z.-W. Sun, Super congruences and Euler numbers, Sci. China Math. 54 (2011) 
2509-2535; preprint larXiv: 1001 .4453^ 19 [math. NT] (2011). 

[22] Z.-W. Sun, On congruences related to central binomial coefficients, J. 
Number Theory 131 (2011) 2219-2238; preprint larXiv : 9 1 1 . 2 4 1 5fr l 6 
[math. NT] (2011). 

[23] Z.-W. Sun, A refinement of a congruence result by van Hamme and Mortenson, 
accepted for publication in Illinois J. Math.; preprint arXiv: 1011 . 1902^ 5 
[math. NT] (2011). 

[24] H.S. Vandiver, Note on Euler number criteria for the first case of Fermat's last 
theorem, Amer. J. Math. 62 (1940) 79-82. 



10 



